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I. INTRODUCTION 

Since its discovery in 1990 J3 the phenomenon of oscil- 
latory interlayer exchange coupling has been the subject 
of an intense experimental and theoretical activity.Ha Its 
mechanism is now well understood and can be attributed 
to the quantum size effect due to spin-dependent confmc- 
ment of electrons in the non-magnetic spacer layer n~U 

A remarkable result of the theory of interlayer ex- 
change coupling is that it becomes particularly simple in 
the limit of large spacer thickness D (asymptotic regime): 
(i) the periods of oscillations versus spacer layer thickness 
are uniquely determined by stationary spanning, vectors 
of the bulk Fermi surface of the spacer material,tj (ii) the 
oscillation amplitude has a universal D~ 2 decayO (ex- 
cept for special nested Fermi, surfaces which give rise to 
a D~ 3 / 2 or D^ 1 decay lawEi), (iii) the amplitudes and 
phases of the oscillatory components are determined re- 
spectively by the modules and arguments of the spin 
asymmetry of the (complex) electron reflection coeffi- 
cients, on the ferromagnetic layers bounding the spacer 
layer .HTD 

Although these simple rules stricly speaking hold only 
in the limit of infinite spacer thickness, they proved to 
be extremely successful in explaining and predicting the 
interlayer exchange coupling observed experimentally for 
spacep-,tiiicknesses down to less than 10 atomic layers 
(AL).0'E3 The surprisingly good results of the asymp- 
totic approximation are also confirmed by first-principles 
calculations .EjTLS 

The exact expression of the interlayer exchange cou- 
pling may be expanded in a series of powers of 1/D. The 
first non-zero term is the D~ 2 term; it corresponds to 
the asymptotic approximation. The higher order terms 
(varying like D~ n with n > 3) are called the preasymp- 
totic corrections. 

The overall success of the asymptotic approximation 
in explaining the experimental results obtained at fairly 
low spacer thicknesses is a priori surprising and raises 
the question of a quantitative determination of the range 
of validity of the asymptotic approximation. In addition, 
there are also some cases where the asymptotic approx- 
imation fails in providing an accurate description of the 
interlayer coupling for small spacer thicknesses. 



The inadequacy of the asymptotic approximation for 
an accurate description of the interlayer exchange cou- 
pling is best illustrated by the Co/Cu/Co(001) sys- 
tem. The "exact" calculations (i.e., not relying on the 
asymptotic approximation) show that the D~ 2 decay 
law of short period oscillation is not obeyed .until the 
spacer thickness reaches approximately 20 AL.E3 In ad- 
dition, the amplitude of the long period oscillation for 
Co/Cu/Co(001) with thick CoJayars as calculated from 
the asymptotic approximationEiffii is typically at least 
one order of magnitude too small as compared to the 
one obtained from calculations . wh ich do not rely on 
the asymptotic approximation.t§tZl The inadequacy of 
the asymptotic approximation to describe accurately the 
long period oscillation of the Co/Cu/Co(001) system 
with thick Co layers is also illustrated by the fact that 
the corresponding amplitude, as well as the period itself 
depends o p - th e Cu thickness range which is used to de- 
termine itE-Da and usually does not agree jsdth the one 
calculated from the bulk Cu Fermi surface. cB 

The present paper is devoted to a detailled discussion 
of the validity of the asymptotic approximation and of 
the preasymptotic corrections. It is organized as fol- 
lows. First, in Sec. || I recall the assumptions under- 
lying the asymptotic approximation, and the results that 
it yields. Next, in Sec. |lj I present a systematic method 
for calculating the preasymptotic corrections to arbi- 
trary order and their temperature dependence, and we 
carry out explicitely the calculation of the first preasymp- 



totic correction, i.e., the D term. The Sec. IV is de- 



voted to the discussion of the results with emphasis on 
the Co/Cu/Co(001) system, and concluding remarks are 
given in Sec. M. 
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II. ASYMPTOTIC APPROXIMATION 

he general form of the interlayer exchange coupling 
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In this expression, g(e,kii) = fc+ — k z is the difference with the complex amplitude determined by 



between wave- vectors of an electron propagating through 
the spacer layer in the +z and —z directions (the z axis 
is taken perpendicular to the layer plane). Here a sin- 
gle contribution has been considered; in the general case 
there would be several such contributions, due to multi- 
ple bands in the spacer materials and to higher harmon- 
ics (i.e., higher order terms in an expansion in powers 
of e 19 ' e ' k "' ) ) but the calculation of the various contribu- 
tions is exactly the same; thus, for the sake of simplicity, 
a single contribution is considered here. As explicitely 
indicated, g(e, k||) varies with the energy e and with the 
in-plane wave- vector || . The other factor in the integrand 
of Eq. ([!]) are the Fermi-Dirac function f(e,T) and the 
complex amplitude M(e, k||) which depends on the spin- 
asymmetry of the reflection coefficients at the spacer- 
ferromagnet interfaces.0 The range of the integration over 
k|| in Eq. (Q) is the surface Brillouin zone corresponding 
to the crystalline orientation of the layers. 

The asymptotic approximation is based upon the ob- 
servation that, because of the strong variation of the 
Fermi-Dirac function at the Fermi energy, and because 
of the rapid variation with e and k|| of the exponential 
factor, the behavior of Eq. ([!]) at large D is dominated 
by the contribution of states on the Fermi surface, such 
that the spanning vector of the Eermi surface, q{e, kn), 
is stationary with respect to k|| a In the general case, 
there may several such stationary spanning vectors, each 
of them giving rise to an oscillatory component of the 
interlayer exchange coupling, but for simplicity, I shall 
consider the case of a single component; the generaliza- 
tion to the case of multiple components is immediate. 

The in-plane wave- vector kjj corresponding to the sta- 
tionary spanning vector of the Fermi surface q* is taken 
as the origin for k|| and the Fermi level is taken as the 
origin for s. 

The variation of M(e, kn) with e and k|| is neglected, 
i.e., we assume 
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The temperature dependence is given by the function 
F (x) 



sinhx 
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The predictions of the asymptotic approximation are 
summarized as follows: (i) at T = the coupling is given 
by a periodic function of D multiplied by a decay factor 
D~ 2 , (ii) the phase of the oscillations is independent of 
the temperature, (iii) the complex amplitude factor varies 
with spacer thickness and temperature only as a function 
of the product DT. 



ILL CALCULATION OF THE PREASYMPTOTIC 
CORRECTIONS 

Let us compute the preasymptotic corrections that 
appear when we release the approximations made in 
the previous section. First, instead of assuming that 
M(e, kn) is a constant, we expand it around e F = and 
kjj ee as 
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which defines the expansion coefficients M p . q ^ r . In addi- 
tion, we extend the expansion of q(e, k|| ) beyond the first 
order in e and the second order in k x and k y : 

q(s, k„) = q* + ^- - - ^ + Q(e, k„ ) (8) 
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where Mo ; o,o is a constant (the motivation for the choice 
of the notations will appear clearly below). Further, we 
expand q(e,\c») near q* as 



g(e,ki 



2e 
hvp 



k x 

K.i 



2v_ 

Kn 



(3) 



where Dp is the Fermi velocity, and k x and k v the curva- 
ture radii of the Fermi surface corresponding to kjj (the 
x and y axes are chosen so as to eliminate the term pro- 
portional to k x k y ). Finally, since only the neighbohood 
of kjj contributes significantly to the integral, we extend 
the integration range for k x and k y from — oo to +oo. 
Then, we obtain easilyll 
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which defines the expansion coefficients Q s -t. u - The fac- 
tor e lqD (the e and ky arguments will be dropped in the 
following) in Eq. (0) is then rewritten as 
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Then the product M(e, k|| with the last factor of the 
above equation is expanded in powers of e, k x and k y as 
follows 
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which defines the new expansion coefficients C p:q ^(D). 
Their explicit expression in terms of the coefficients 
Mp. qjr and Q S ;t,u can be obtained by a straighforward 
term-by-term identification. 

Inserting the above expressions in Eq. ([!]), we obtain 
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The above expression is non-zero only if q and r are 
even. Performing the integrations as explained in the 
Appendix, we obtain 
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and where the functions F n (x) governing the temperature 
dependence of the coupling are defined by 



dx n V sinh x 
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Finally, we reorder the terms of Eq. ( |l3| ) so as to express 
it as a series in powers of 1/D: 
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The expressions of the functions A n (DT) in terms of the 
coefficients M p - qjr and Q s -,t,u are obtained from Eq. Jl3| ) 
and from the explicit expression of the expansion coef- 
ficients Cp. q ^ r (D) after tedious but straightforward alge- 
braic manipulations. Below, I give the explicit expression 
of the asymptotic approximation term Aq and of the first 
preasymptotic correction term A\ : 
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IV. DISCUSSION 
A. Symmetry considerations 

As discussed in detail in Ref. |?], the various critical 
points k| of the two-dimensional Brillouin zone giving 
the oscillatory contributions to the interlayer exchange 
coupling can classified according to their symmetry: fol- 
lowing the terminology introduced in Ref. (?], critical 
points corresponding to high-symmetry points of the two- 
dimensional Brillouin zone are termed essential, critical 
points located on high-symmetry lines are termed semi- 
essential, while critical points possessing no particular 
symmetry are termed accidental. 

The symmetry considerations impose a number of re- 
strictions on the c oefficients M p . }q>r , Q s -,t,u and C p - q>r . 
These are given below: 

(i) For an essential critical point, we have M p . qir = 0, 
Qp,q,r = and C P;9 ,r = if q or r is odd. If in addi- 
tion, the critical point corresponds to a rotation axis of 
order equal to or larger than 3, then k x = K y , M p; 2 q ,2r = 

M p -2r,2qi Qp:2q,2r = Qp:2r,2q, and C p -2q,2r = Cp;2r,2g- 

(ii) For a semi-essential critical point (without restric- 
tion, we chose the y axis parallel to the high-symmetry 
line), we have M p . qr = 0, Q p qr — and C p - qr = if 
q = 0. 

(iii) For an accidental critical point, no restriction ap- 
plies. 



B. Preasymptotic corrections at T — 

If we retain only the first preasymptotic correction 
term in Eq. (^), the coupling may be expressed as 
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where the complex length 
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characterizes the preasymptotic correction. Thus, for 
D > |A|, one has 
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where <p and 8 are respectively the arguments of Aq and 
A. From the above equation, it appears that the real part 
of A yields a correction of the amplitude of the oscillatory 
coupling; depending on whether cos 8 is positive or nega- 
tive, one has an increase or a decrease of the amplitude. 
On the other hand, the imaginary part of A contributes 
to a correction of the phase of the oscillation; i.e., for low 
spacer thickness, the apparent period of the oscillations 
differs from the asymptotic onec3 



and becomes 
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depending on whether sin<5 and q* are of the same sign 
or not, one as a decrease or an increase of the oscillation 
period. 

Finally, the above discussion shows that a quantitative 
criterion for the validity of the asymptotic approximation 
is given by two following conditions 
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C. Temperature dependence of the preasymptotic 
corrections 

The temperature dependence of the coupling is gov- 
erned by the functions F n (x) defined by Eq. (J15|). The 
functions F n (x) with n = - 3 are given explicitely be- 
low: 
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FIG. 1. Plot of the functions F n {x) with n = - 3. 

They are displayed on Fig. [y. The functions F n (x) have 
the following general properties: 



F n (x) Mil I3 n x 2 for x < 1, 
F n (x) w - 



x n+1 e x for x > n + 1, 



where the constants f3 n are given by 
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F n (x) with n > 1 has a maximum at 
x* n + 1 , and 
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Because of the marked difference between Fq(x) and 
the functions F n (x) with n > 1, the temperature depen- 
dence of the preasymptotic corrections is in general not 
the same as the one of the asymptotic term. In other 
words, the parameter A char acter ized the preasymptotic 
corrections defined in Sec. [VE generally varies with 
the temperature; thus, the extension of the asymptotic 
regime may depend on the t emp erature. As expected 
and appears clearly from Eq. (17t), the temperature de- 
pendence of A arises from the terms proportional to _Fi 
and i<2, i-e. from the corrections due to the e variation 
of M(e,ku) and Q(e, kii). 



D. Comparison with previous work 

The question of preasymptotic corrections-has been ad- 
dressed previously by Mathon et al£3lE3'c3 They have 
considered the correction due to the energy dependence 
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of the argument of M(e, k||). This correspond to assume 
that the only non-zero correction parameter is -Mi ; o,o and 
that the ratio Mi ; o,o/-Mo ; o,o is purely imaginary. In this 
case, they find that the preasymptotic corrections mod- 
ify the amplitude and the temperature dependence of the 
coupling. This conclusion is in agreement with the result 
of the present paper, if we restrict our analysis to the 
same assumptions. However, even within the restrictive 
assumptions, the method used by Mathon et al. is cor- 
rect only to the order D~ 3 because the the terms order 
D~ n with (n > 4) which are neglected are of the same 
magnitude as the ones which are included. 

Furthermore, as the present analysis shows, there are 
other sources of preasymptotic corrections, for example 
the k|| variation of M(e,k||), that a priori, are equally 
important as the ones considered by Mathon et al. 



E. Discussion of the Co/Cu/Co(001) case 

We shall now focus on the particular case of the 
Co/Cu/Co(001) system, which is often considered as a 
model system for the problem of interlayer exchange cou- 
pling. 



1. Summary of the available data 

Let us first recall the predictions of the asymptotic ap- 
proximation for this system. The most striking result of 
the asymptotic approximation here is that the interlayer 
exchange coupling comprises a long-period oscillation 
originating for the center T of the two-dimensional Bril- 
louin zone (essential critical point) and of a short-period 
oscillation originating from 4 equivalent, semi-essential 
critical points located on the T-X linespQ this prediction 
is confirmed unambiguously both by experimentsEH and 
by "exact" (i.e., not mJviHg on the asymptotic approx- 
imation) calculations,E5~liZl and the periods are in good 
agreement with the ones predicted by the asymptotic ap- 
proximation. In the following, quantities related to the 
short-period and long-period components will labelled by 
S and L indices, respectively. 

Quite generally, the function M(s^\^) determining the 
strength of the coupling is given by 
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where r' (r^) is the reflection coefficient for electrons 
with the spin parallel (antiparallel) to the majority spin 
of Co. Since the majority spin band structure of Co is 
very close to the Cu band structure, one has r^ « and 
it is sufficient to consider r- 1 . 

For the short period oscillation, one has total reflec- 
tion for minority spin electron, due to a local gap (for the 
bands of relevant symmetry) in the minority spin band 



structure of Co, i.e., |r^| = 1.II3~E3 As a consequence, 
the amplitude of the short period oscillation is large and 
depends only very weakly upon the C thickness. These 
facts are confirmed by "exact" calculations .E^TL] The pe- 
riod of oscillation obtained from an "exact" calculation, 
(As = 2.53 AL)EJ is in excellent agreement with the one 
calculated within the asymptotic approximation with the 
same set of potential parameters (As — 2.50 AL)j£j fur- 
thermore, the apparent period agrees very well with the 
asymptotic one down to Cu thicknesses as low a 5 AL. On 
the other hand, as appears very clearly from the Fig. 3a 
of Ref. [l7| the asymptotic D~ 2 decay law is satisfied 
only for Cu spacer thicknesses larger 20 AL; in practice, 
this means that for a quantitative comparison with ex- 
periment results (which are usually obtained for spacer 
ticknesses smaller than 20 AL), the amplitude calculated 
from the asymptotic approximation is inappropriate. 

For the long period oscillation, the situation is quite 
different. Here, the reflection coefficient rl is miite small 
for a semi-infinite Co layer (|r£| » 0.15)ill@B@ which 
leads to a very small amplitude for the long-period os- 
cillatory component; all theoreticaJ calculation, whether 
"asymptotic"ll5L3 or "exact" E£tl3 agree on this point. 
In addition, there is a strong variation of the amplitude 
of the long period oscillation with Co layer thickness, 
which is due to quantum interferences in the Co layer, — 
a prediction which has been confirmed ^pe,rimentallyL2 
and by "exact" theoreticaL-Gak«lations.ExE-3^ However, 
although both asymptotidl3~E3 and "exact"liinl3 calcu- 
lations agree on the fact that the long period oscilla- 
tion has a very weak amplitude for thick Co layers, 
they disagree on the value of the amplitude: while "ex- 
act" calculation yield an amplitude of the long period 
oscillation -of—the order of 10 % of the short period 
amplitudeJiSEJjO asymptotic calculations yield a ratio of 
the order of 1 % onlyO~E3 In addition, the value of the 
long period coupling obtained from "exact" calculations, 
Al = 5.09 AL,El differs markedly from the one calculated 
from the asymptotic approximation with the same set of 
potential parameters, Al = 6.49 ALpi and depends on 
the thickness range which is used to determine it as is 
seen from Fig. 2 of Ref. [IT]. 

So, our purpose is to explain why the asymptotic ap- 
proximation works well for some aspects of the coupling 
and departs markedly from "exact" calculations for some 
other aspects. 

For both the long-period and short-period contribu- 
tion, the preasymptotic corrections due to the e and k|| 
dependence of Q(e, kn) is unimportant, so that we shall 
focus on the corrections due the variation of M(e, k||). 
Thus, by taking the symmetry considerations mentioned 
in Sec. IV A, the corrections are given by 
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respectively. Let us now recall that within elementary 
models such as the free-electron model or the single-band 
tight-binding model, the reflection coefficient for a semi- 
infine barrier is either a real number of module smaller 
than f (partial reflection) or a complex number of mod- 
ule equal to 1 (total reflection). As shown in Ref. this 
remains approximately true for realistic multi-band sys- 
tems under rather general circumstances. 



2. Short-period oscillation 

Thus, for the short period oscillation, only the phase 
of the reflection coefficient varies with e and k|| , and 
as a consequence, As is a real n umbe r. Hence , ac- 
cording to the discussion of Sees. IV B| and IV C , the 
preasymptotic correction affects only the amplitude of 
the coupling and its temperature dependence, while the 
period is not affected by the preasymptotic corrections; 
these conclusions are in agreement with the results sum- 
marized above. Similar conclusions have been obtained 
by Mathon et al!E3 who have considered the preasymp- 
totic correction associated with the energy dependence 
Af(e,k||), i.e., with Mq. q1 and neglected the one associ- 
ated with its k|| dependence, i.e., with Mq. 2 and Mq. 2 ; 
indeed in view of the similarity of their corrected results 
(see Fig. 13 of Ref. E3) with those of "exact" calculations 
(see Fig. 3 of Ref. it seems that the Mq. x term al- 
ready accounts for a large part of the total preasymptotic 
correction. Further work would be needed to assess the 
importance of the Mq. 2 and Mq. 2 terms. 



than the asymptotic one; this conclusion provides an con- 
sistent explanation for the discrepancy on the long oscil- 
lation period mentioned above. However, the preasymp- 
totic regime here has an unusually large extension and 
asymptotic behavior is expecte d to h old only for D ^> 
50 AL; thus, the analysis of Sec. IV B , which was limited 
to the lowest order in 1/D, is certainly not sufficient to 
analyse the results obtained from experiments and from 
"exact" calculations and it is not surprising that not only 
the period itself, but also the amplitude of the long period 
oscillation predicted from the asymptotic approximation 
is inappropriate for spacer thicknesses D < 50 AL. 

Actually, 90% of the result @ for A L is due to the ky 
dependence of r^, so that the unusually long preasymp- 
totic regime for the long-period oscillation is to be at- 
tributed to the presence of a narrow window of low re- 
flectivity near T in a region of otherwise total reflectiv- 
ity, for minority spin. This pecularity can in turn be 
explained by a rather simple following argument: The 
relevant band in the Cu is mostly of p z character. The 
reflectivity is due to the hybridization of the correspond- 
ing p z band in Co with the minority spin d bands. At 
the center T of the two-dimensional Brillouin zone, only 
the d^ z 2_ r 2 band is allow by symmetry to hybridize with 
the p z band, and hence to contribute to the reflectivity; 
however, because the d\ z i_ r i band is almost full, it yields 
only a weak reflection coefficient at the Fermi level, and 
total reflection is attained only 0.5 eV below ep- But, 
as one moves away for T, the d\. z and d^ z bands that 
lie close to the Fermi level are allowed to hybridize with 
the p z band, which yields a strong increase of and, 
eventually, total reflection. 



V. CONCLUSION 



3. Long-period oscillation 



For the long-period oscillation, on the other hand, since 
r L is real, only the magnitude varies with e and k|| and 
the phase is constant; thus, At, is purely imaginary. As 
discussed above, this leads to an apparent shift of the pe- 
riod in the preasymptotic range, which is precisely what 
happens in this case. Let us attempt to estimate the 
value of At,. As seen from, e.g., Fig. 22 of Ref. ^ r± 
increases with decreasing energy. Furthermore, as seen 



from Fig. 2 of Ref. 18 and Fig. 2 of Ref. El, r* increases 



very strongly with k|| and full reflection is reached at a 
distance 0.1 x tt/cl from T; indeed, the low reflectivity 
arises only in a narrow window around T. Taking these 
two contributions into account, we arrive at the result: 



A, 



-i x 50 AL. 



(30) 



Because q* L is negative, this result implies that the 
preasymptotic correction will lead to an apparent oscil- 
lation period in the preasymptotic regime that is shorter 



I have presented a detailed discussion of the question 
of preasymptotic corrections for the intcrlaycr exchange 
coupling. A systematic method for computing exactly 
the preasymptotic corrections to arbitrary order in D^ 1 
has been presented, and the explicit expression of the 
first correction term (of order D~ 3 ) has been given. 

This method allows one to assess quantitatively the 
spacer thickness range in which the asymptotic approxi- 
mation is expected to be reliable, and the one in which 
preasymptotic corrections should be taken into account. 
In the latter case, I have shown that the preasymptotic 
correction alters the amplitude and/or the apparent os- 
cillation period, depending on the argument of the cor- 
rection parameter A. 

The case of Co/Cu/Co(001) has been discussed in de- 
tail. I have shown that most of the discrepancy between 
the asymptotic approximation and "exact" calculations 
for the system can be understood on the basis of the the- 
ory for preasymptotic corrections presented in this paper. 
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APPENDIX: 
Integration over energy 

We consider here the following integral 

" + OO 



/-t-oc 
de/(e,T)e n exp (2iea) . (Al) 
-oo 

We can compute it by using the method of residues. The 
Fermi-Dirac function f(e,T) has poles for 



s p = m{2p + 1) 0- 



(A2) 



where p is an integer and (3 = (fc^T) -1 . The correspond- 
ing residues are equal to — /3 _1 . Closing the integration 
path in the upper half of the complex plane, we obtain 

I n (T) = -2[j\ ^(2p+l) n e -^p+^/P. (A3) 



Finally, we obtain 



p=0 



I n (T) = 7° F n 



2n a 



whith 



and 



7° = -i n+1 a- (n+1 ^n\ 



= (-l) ra x n +1 d" ( I 
n\ Ax n Vsinha; 



/•',.(.'•} = - — 7— ' 



(A4) 



(A5) 



(A6) 



Integration over wave vector 

Here we calculate the following integral. 

/■oo 

B 2n = dk k 2n exp (-i7fc 2 ) . (A7) 



We can deform the integration path in the complex plane 
in such a way that the point k ^-,0 be traversed along 
the direction of steepest descent .Eil The integral is then 
easily calculated, and we get: 



Bo 



1/2 (2n + l)!! 
2™ (2n+ 1) 



(A8) 
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